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DIFFUSION AS A FUNCTION OF AGGREGATION IN 
COLLOIDAL MEDIA 


HERMAN BRANSON 
HOWARD UNIVERSITY, WASHINGTON, D. C. 


With the assumption that adsorption is a simple function of surface 
area, an analytical treatment is given for the dependence of the diffusion 
coefficient of an adsorbable solute upon the degree of aggregation of the 
adsorbing colloid. A simple relation is deduced after introducing some 
approximations. Some implications of the final diffusion equation are 
given. 


The effects which colloids may have upon diffusion have been 
discussed qualitatively by N. Rashevsky (1938) and the specific prob- 
lem of the adsorption of a diffusing solute by a single type of colloid 
has been treated by J. M. Reiner (1939). This note reports an investi- 
gation of a general phenomenon in a colloidal medium: with the 
assumption that adsorption is a simple function of surface area what 
could we expect to be the effect upon the diffusion coefficient of an 
adsorbable solute of the continued aggregation of the adsorbing 
colloid. This discussion proceeds under some plausible assumptions 
which simplify the analysis. The results are finally related to certain 
membrane phenomena of the cell. 

The system we consider has initially simple, colloid particles and 
solute particles. The primordial colloid particles coalesce giving larger 
particles of uniform surface, uniformly distributed. If the initial 
number of simple colloid particles is designated by N- and the initial 
number of neutral solute by N,, the relationships which must always 


be met are 


N,=2, + 3S In(ha,k) 
=i) Re2 


(1) 


M2 


N.= 


S k[n(c, k) + n(1,4,k)] , 


l=1 


where n(l, a, k) designates the number of k-th-aggregated colloid par- 
ticles with J adsorbed solute particles, n(c, k) is the k-th-aggregated 
colloid particles having no adsorbed solute particles, and ”, is the 
number of solute particles not adsorbed. The types of possible re- 
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actions taking place in the system may be written 
n(c,1) +n, > n(1, a, 1) 
n(1, a, 1) +(e; 1) > We are) 
Mee ray ths eS (2) 

n(c,i) +n(c,j7) >n(ce,i+ 7) 

n(i,a;7) +n, 2 ni +1, 4,7) 


provided i+ 17, This condition is explained below. The first 
three reactions express the general assumption that a solute particle 
may be adsorbed by a single colloid particle but that the aggregate 
cannot adsorb another solute particle until the surface is at least 
twice the surface of the single colloid particle. Any other reaction of 
agglomeration is possible in the system so long as it does not violate 
this principle. 

Nevertheless two particles having such a number of adsorbed 
particles that the complex formed by their agglomeration would have 
more than the allowed number of adsorbed solute particles were they 
all retained, can still agglomerate with the release of solute particles. 
This is expressed by adding to (2) 

WAT, OL) eC Os Le Siete 


n(k,a,m) +n(j,a,7r) >~n(q,a,m+r7r) + In, 


where (m+71)?82—1<qS(m+r)*, ktjg=qt+l. 

For simplicity the reactions are assumed to be irreversible. With 
these assumptions the maximum number of adsorbed particles is 
easily deduced. If we let v,, 7,, and s; represent the volume, radius 
and surface of a single colloid particle and v;, 7, and s; the same for 
an aggregate of k particles, we have 


4 
i= 3 te > S:—4AnrZ 


VE = kv, 


Since the k aggregate to form a spherical cluster, kr,2—= 7,2 thereupon 
Si Ank?/*r,? = k?/3s,, The maximum number of adsorbed solute par- 
ticles on an aggregate of & colloid particles is 


l=)" if k is an integer, 


ee <p aes if k is not an integer, 
which we have already used in (2) and (8). 
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Some hint as to possible numerical values for the maximum num- 
ber of colloidal particles which may aggregate before flocculation can 
be taken from the specification of the range of colloidal particle size, 
the diameters ranging from 10~ to 10° em. Assuming that our prim- 
oridal colloid particle has the smaller diameter and the largest the 
larger, we have the maximum for mo 10° and the maximum for 
Ecol 0+: 

The solute particles are taken as uncharged so Fick’s law for the 
transport of solute is 


J=JI,+-3 SI (1,0, k), (4) 


in which J is the number of solute particles passing through unit sur- 
face in unit time whether as simple solute or in aggregate. The/, k, 
r,and m are the same as in (1). The J(1,a,k) may be written 


J(l,a,k) =— D(l,a,k) gradn(l,a,k). (5) 


Introducing (5) into (4) with the corresponding forms for J; and 
J(l,a,k) and with the dependence of n, and n(l, a, &) on the posi- 
tional coordinates expressed through N, , we have 


— D grad N; —— D, grad 7,:— > 3 Dil, a,k) grad n(l, a, k) 


t=1 k=1 
or 


ON: ae a on(l,@,k) 
Sara eA Aol) saan ae an 


13 — 3 
ON, t=1 kea1. ONs 


(6) 


D is the effective diffusion coefficient of the solute. Using (1) this be- 
comes 


Sete (1 a Je) 7m on(l,a,k) 
= — | ——_—— |+ Dil, a, k) ——————_ (7) 
D=D,|1— 3 BI #2 22G 2m Sy 


Inasmuch as the adsorbed solute particles do not appreciably affect 
the bulk of the colloid aggregates we may take D(l, a, k) = D(c, k) 
for all values of k. Furthermore on the basis of Einstein’s law of dif- 
fusion the diffusion coefficient of an aggregate of h colloid particles is 
related to that of an aggregate of one solute particle and one solute 
particle by 

D(e ,1) 


Di(c,h) rake hala Tis 2 (8) 
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These substitutions simplify (6). There remains, however, the rather 
complex double summations. We can approximate them through 


¥ Sin(1, a, k) © B(p)n(7, 4, p) (9) 
SE Sn(l,a, kb) oa(p)n(r, 4, p). (10) 


Without reference to the coefficients or to the first and last symbols 
in the parentheses we see that the n(l, a, k) are merely the number 
of particles of each size group in the system. The mechanism here 
described places all the particles initially in n(1, 4,1), thus initially 
(10) could be replaced by n(1, a, 1). As time goes on, other 
n(l, a, k)’s will arise but the most unlikely distribution, one per- 
fectly flat, would lead to a multiplicative factor equal to the number 
of particle sizes. For any other distribution the factor for any p would 
be in general a function of p but less than p. The factor 6(p) would 
behave in the same manner, although the numerical value would be 
different from that of a(y). We can relieve ourselves from too much 
specificity by putting 


a(p) =p",a<1;8(p) =p’, b=1. 


On closer inspection, (10) is seen to be equivalent to replacing the 
size-distribution curve by a rectangle of height n(7, a, ») and width 
a(p). In (9) we replace the three dimensional distribution by a solid 
of cross-sectional area $(p) and height n(r, a, p). It is of interest 
to observe that » will be a function of time, for the distribution curve 
of particle size will move towards the maximum size with increasing 
time. Instead of the variation of n(r, a, ») with time we think it 
best to keep p, which we shall call the degree of aggregation, as the 
independent variable. 

If N, > N. then practically all the N, will be used as n(1, a, 1) 
before the formation of any higher associated particles. This permits 
the division of the problem into two sections, the initial aggregation 
to form n(1, a, 1) and the subsequent reactions to form higher ag- 
gregates. This first is essentially J. M. Reiner’s problem. Both are 


contained in our analysis as we see upon introducing the approxima- 
tions and the relations 


n(1,a,1) N- 
n(l,a,p) © ———_ = — 
p 


and —N,=@N, vere 
whence 
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D=D,[1— pC?) D(1,a,1) dn(l,a,p) 
oN, (2/5 0 ) oN, 
Rte | € peel) (1s, 05/1) Cy 
‘| aoeN 25/3 —b) ; 


In terms of time, (11) holds from ¢, to tn where (t, — t) would be 
the time required for all N. to become n(1,a,1) and t, — t; would 
be the time for the formation of maximum size. For t, S t S t, Rein- 
er’s (3) holds 


On Cl = By onl, at) 
D=—D, a a ie ee eee 
oN, D(1,a,1) oN, ? 


7) 
PERMEABILITY 
h— 


DIFFUSION COEFFICIENT 9 
oD—~ 


o 


11,0,0) 


i 
DEGREE OF AGGREGATION 
Pp 


FIGURE 1 


which is observed to be the same as our (11) with p = 1. Repre- 
sented graphically these equations join exactly at p = 1 as is shown 
in Figure 1. 

The removal of the assumption that N, > N, and the precise se- 
quential reactions do not influence the shape of the curve. This would 
not contradict (3) so long as the binding energy of colloid for colloid 
is not considerably greater than that of colloid for solute. It may 
cause a smoother transition at the minimum and a raising of the mini- 
mum so that D does not reach D(1,a@,1). Furthermore we need not 
limit the number of adsorbed particles to one for the surface area of 
a simple colloid particle. If we would consider the number of adsorbed 
particles per cm?, then the resulting curve for the variation of the dif- 
fusion coefficient would probably have the same shape. We may also 
weaken the requirement that the agglomeration gives spherical par- 
ticles by introducing a shape factor y, so that 


Ss=yhts,, bs >yZ1. 


The shape factor will be a function of k, the number of particles in 
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the complex. With our assumptions the inclusion of this factor would 
raise the maximum for 1 with a given k to y k#?? —1<1S y k** but 
would cause no change in the other equations. 

E. L. MacBain’s data on the diffusion coefficient of long chain 
sulfonic acids, C1. to C,,, exhibit a minimum at dilute solution fol- 
lowed by a rise at higher concentrations. The curve is similar to 
Figure 1 but the data do not seem susceptible to interpretation by our 
mechanism These experiments are instructive, however, in calling 
our attention to the fact that the D-p curve in Figure 1 could be trans- 
ferred into a D-concentration curve. Since the aggregation is directly 
proportional to the concentrations by measuring the diffusion coeffici- 
ent of various concentrations of N, and N, a test of the mechanism 
outlined here could be carried through so long as the degree of ag- 
gregation does not vary appreciably during an experimental deter- 
mination. 

If the reaction takes place within a cell we may take N, = q,¢ and 
N. = qcet where the q’s are the average rates of production of the so- 
lute and the adsorbing particle within the cell. In (11) e« would be 
Q-/ds. For sufficiently rapid aggregation the action of this mechan- 
ism within a cell may give a purely physical interpretation of certain 
trigger mechanisms and other aspects of cellular behavior. 

The rather simple extension that the presence of ”, primes the 
release of n, in spurts gives a sinusoidal curve for the variation of 
the diffusion coefficient. 

Finally the view that the diffusion solute lowers the permeability 
of the cell membrane to another diffusing substance admits a varia- 
tion of the membrane permeability, h, inverse to that of D in Fig- 
ure 1. 

The author wishes to express his thanks to the Julius Rosenwald 
Fund for a grant which supported this investigation. He wishes also 
to acknowledge the helpful suggestions of Prof. N. Rashevsky out of 
which came this problem and the discussions with the University of 
Chicago group of mathematical biophysicists. The author is especial- 
ly indebted to Dr. Alvin Weinberg who read the paper and made some 
valuable suggestions. 
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A THEORY OF STEADY-STATE ACTIVITY IN NERVE-FIBER 
NETWORKS: IV. N CIRCUITS WITH A COMMON SYNAPSE 


ALSTON S. HOUSEHOLDER 
THE UNIVERSITY OF CHICAGO 


Conditions under which either of two distinct activity patterns may 
arise from the same stimulus pattern are deduced for the case of a net- 
work which consists of N simple circuits all joined at a common synapse. 
If the product of the activity parameters of all the fibers in any circuit 
is called the activity parameter of the circuit, or, more briefly, the cir- 
cuit parameter, then the condition for the existence of such mutually 
consistent activity patterns is that there be a sum of circuit parameters 
which is not less than unity. 


The general problem which is being discussed is twofold: first, to 
characterize those fiber networks, in which the activity pattern (AP) 
is uniquely determined by the stimulus pattern (SP) ; and second, for 
those networks which do not. possess this unicity property, to char- 
acterize those pairs of AP’s which are consistent with the same SP. 
In previous papers (Householder, 1941; these papers will be referred 
to here as I, II and III) this problem was solved for the simple-circuit 
case, and reference is to be made to them for these results and for the 
assumptions and terminology here employed. We remark only that 
the linear case alone (steady-state response a linear function of the 
constant stimulus) is being considered at present. 

In working toward the case of the most general possible net- 
works, we proceed now to the case of next greater generality beyond 
that of the simple circuit. We consider a network which consists of 
N circuits. The fibers in the J-th circuit are represented by 


(1), (23), +++, (ms). (1) 
Here and hereafter, when not further restricted, the indices 
yd ply ee ye 5 IN (2) 
and 
tgs Jrytet HAs, MW. (3) 


The sub-index J will be omitted in cases where it is evident from the 
context, or where it is not necessary to make the specification. All 
fibers n,; terminate together and synapse with all fibers (1x), thus 
connecting all the circuits at this synapse. Moreover in any circuit 


7 


8 NERVE FIBER NETWORKS 


J , the first fiber, the first two fibers, --- , or the first n; — 1 fibers may 
be identical, fiber for fiber, with the corresponding first fiber or fibers 


4 ( ) 22 
as 
2. 3 
Sy 
FIGURE 1 


of some other circuit, but beyond these the two circuits have neither 
fibers nor synapses in common. The figure pictures such a set of three 
circuits in which the fibers (1.) and (1;) are identical. 

In II we classified all fibers of a circuit into “even” and “odd”. 
This classification can be made as follows: the first fiber in the cir- 
cuit is odd; the i,-th fiber is odd provided the product of the activity 
parameters of the fibers which precede it in the circwt is positive; 
otherwise the fiber is even. In the present case, the “first” fiber is 
determined by the structure of the system, and the classification is 
unambiguous. We shall call an AP “even” (“odd’) for a given circuit 
or set of circuits provided all its zeroes in that circuit, or those cir- 
cuits are at even (odd) fibers. 

The quantities y;, z:, o; can be defined exactly as before, with 
the understanding that some of these quantities are identical when 
there are fibers common to two or more circuits, and that 


Si, ot Si, (4) 


for every J and K. Also it is, of course, tacitly understood throughout 
that the possible AP’s are always limited by the fact that the same 
stimulus acts upon all the “first” fibers, so that when, for example, 
fiber (1,) is active, every fiber (1x) which has a lower threshold, is 
necessarily also active. 

There is also a further limitation upon the possible AP’s for the 
system. It was shown in III that for a simple circuit it is possible to 
have complete activity only when the product of the activity parame- 
ters of the circuit is less than one (algebraically. Thus complete ac- 
tivity is always possible in a circuit having an odd number of inhibi- 
tory fibers, for the product is negative, and hence less than one). An 
analogous argument can be used to show that 

LEMMA 1. In a multiple circuit of the type under discussion, 
complete activity in the steady state is possible if and only if the prod- 
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uct of the activity parameters of any circuit and every sum of such 
products is less than unity. 
The proof depends upon the fact that if 


Ase cece (5) 
and if 
A;=1-A,, Ayx =1—A,—Ax,-*-, 


A,,...,.vw=1-—D> Ay (6) 


then A,,..., 7 is the determinant of the equations in y; which we obtain 
when we set every a; = @;. 

We therefore confine our attention to AP’s of the following type, 
which we shall call admissible: 

a) if a fiber (1,) is inactive, then every fiber (1x) with the 
same or a higher threshold is also inactive; 

b) if any component simple circuit is completely active, say the 
J-th, then A, < 1, and every sum of such A,’s is less than one. 

We now state the principal result, after which this will be proved 
by a sequence of lemmas. 

THEOREM. In the multiple circuits described above, let a and B 
be distinct admissible AP’s, and let aw, represent the set of zeroes of 
a inthe J-th circuit, Bo, the set of zeroes of B in the J-th circuit. Let the 
circuits be so enumerated that for J =M, a, and fo, have no members 
in common; forM < J =M’, aw, and f., are null sets; for M’< JEN, 
a, and Bo, have at least one common member for each J. Then for the 
existence of ao consistent with both a and f , it is necessary and suf- 
ficient that 

i) every member of a, which is not also a member of fy is even, 
and every member of 8, which is not also a member of a ts odd, or else 
the reverse is true; 

ii) for every J > M’, if j; is a member of both a and fo , and if 
there is an i; which is a member of one but not of both a. and By , then 
i535 

iii) forJ =M, there is at least one A; > 0, and the sum of all 
the positive A; added to Ayu +--- + Am, is not less than one. 

Consider first the case that « and # have no common zero. For 
the sake of more compact expression we introduce the following nota- 


tions: 


for h ea read “his a zero of a”; 


forhea, + 6, read “h is either a zero of a or a zero of £”. 
Let 
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y, when he, 
Wy = (7) 
zr When hea. 
Let 
+1 when h,keo, or h,kefp, 
[hkKj—=[k, h] = (8) 
—1 when hea and Kkef,. 


The w and the bracket symbols are not defined for other values of the 
indices. Let 


A, = Ay Ani ° °° Uy when l=aA< ks Ny. (9) 
LEMMA 2, Ifh, kea),+ fp, and1=h<k=n,;, then 
We = [h, k] Ak wr, (10) 


the symbols being defined by (7), (8) and (9). Except for the sym- 
bols, this is identical with lemma 3 in II, proved there for the simple 
circuit. The same proof applies here. 

In each of the first M circuits there is a member of a + 6.. Let 
h, represent the first of these, k, the last, in the J-th circuit. The case 
h, = k, is not excluded. To further simplify the notation, let 


Yu=Ynx > Zk =2nx » 
Ya — Vix» OR ies 
1k =Vix» Cx = Zkx » (11) 
ia Une, Ox =2nx » 
Wr = Wk, » W k= Wri 


Furthermore, since the discussion is somewhat different for the first 
M and the last N — M circuits, let the indices have the ranges 
BA) Ea gp mei teal" i Hi l= Mila Ne (12) 
while 
L142 ,--*,N. (13) 


Now each fiber (1,) is acted upon by all the active fibers (n,) as 
well as by the external stimulus S,. Hence 


Yi — 01,7 > An, Mea SS Ong Yas 
2, =01, +E Bo, Zy + Edn, Zn, ioe 
so that 
Yt — 2,=—W (15) 
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is independent of L. Also by lemma 1, 

77 — C= Ae (yy — 23) = ele W. (16) 
Since the sign of the left member of (16) depends only upon whether 
h; € o or h,; « 6, this equation together with lemma 2 establishes the 
necessity of condition (i) in the theorem. 


Now by the lemma proved in III, if we subtract the second of 
equations (14) from the first we have 
WD ARH Ue + DAn(Yn— 2x), (17) 
where 
ne i Kee po, 
Ux = (18) 
—f, lf Krea. 


But by (15) and (16), (17) is equivalent to 
Aua,..,.»W=2 As UK » (19) 


and by lemma 1 this A > 0 since a and £ are both admissible. Again 
by the lemma from III, 


ace bp AN ig ages 


so that 
Aged te, w(y’s —C)= Ae > Aveo UK. (20) 


Since the A is positive, each member of (20) is positive when h, «€ Bo 
and otherwise negative. If we change signs in each equation (20) in 
which the left member is negative, and then drop each non-positive 7 
and ¢ from the left (one from each equation), we are left with the in- 
equalities 


Aya,...,yWy AY > (As, kx] Ales Wr. (21) 
Finally we employ lemma 2, which implies in the present case that 
wy, = [hy, ky] Ay W's, (22) 
and combine these inequalities with (21) to obtain 
Avs, --, Nw Wy = B! SL , K] Bx Wk 5 (23) 
where 
Bi Aw, be Ate [dK] [h,, bx) (24) 


(Note that the superscript J is not intended as an exponent.) 

If (21) and (22) are regarded as inequalities for the determina- 
tion of the w’s, then (23) are the necessary and sufficient conditions 
for the consistency of these (cf. Dines, 1919). Likewise, by an argu- 
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ment which parallels that given in II, we can show that (22) and 
(23) are the necessary and sufficient conditions for the determination 
of the remaining w’s in the multiple circuit. If these can be solved, 
then all remaining y’s and z’s, and hence all o’s, are obtainable as so- 
lutions of a consistent system of equations. Hence to show the neces- 
sity and the sufficiency of the conditions of the theorem, we have only 
to show the consistency of (23). We first suppose (23) consistent 
and show the necessity of the conditions of the theorem. 

Now all the w’s are positive, as is the 4. The inequalities (23) 
can also be written in the form 


Aye ce Wy = [1 ’ J] Ba [1 ; K] Bx Wk (25) 


where the sum is independent of J. Hence either [1, J]B’ is positive 
for all J or negative for all J, according as the sum is positive or 
negative. Hence the coefficient of each wx on the right in (25) or in 
(23) has the same sign in every inequality, and is the same as the 
sign of Ax, since this is its coefficient in the K-th inequality. Cer- 
tainly there must be some K for which Ax > 0 if (23) are consistent. 
If, for some K , Ax rss, ..., y = 0, then we have completed our first task 
and shown the necessity of the conditions in question. If not, then we 
can suppose that Ay > 0, and since by assumption Aaas,..., 7 > 0, it 
is legitimate to multiply the J-th inequality (J < M) by Aman,...,», 
the M-th by [J, M]B’By and add. After suppressing the factor 
Axs1,...,~ Which appears on both sides of the result, and is known to 
be positive, there results the system, 


M-1 
Axwu, +, Nw Wy S Be > LJ, K] Bx Wk (J < M) ’ (26) 
K=1 


of M — 1 inequalities independent of wy . 

The argument can be paraphrased for this system with M — 1 
replacing M , and if every Ax,u,...,7 > 0 we obtain a similar system of 
M — 1 inequalities. Proceeding thus sequentially we either reach a 
stage at which there is a negative 4, in which case the necessity of 


the conditions is proved, or else we continue to obtain finally the in- 
equality 


Ags, «.,n Uy S Ay Wy (27) 


as the condition for the consistency of the preceding systems, and this 
can be satisfied only if 


Oe alo) 


This proves the necessity of (iii) in the theorem. 
To prove the sufficiency, if there are negative A,’s, let these be 
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the first M,. For M, < J = M, by hypothesis if we subtract the 
sum of all of the A; from Ayuu,..., 7, the result is negative, while the 
A itself is positive. Hence there must be some maximal number (pos- 
sibly zero) of these A,’s, each taken once, whose sum does not exceed 
this 4. Let this number be M — M, 2 0, and, if M. < M, let Ay, 
Am, -:: be a set of the A’s whose sum does not exceed the 4. We 
may then proceed with our elimination from (23) until we reach a 
set which has the same form as (23) with M, replacing M. In these 
inequalities, when all terms are transposed to the right, the first M, 
of the w’s have everywhere negative coefficients, the remaining 
M, — M, of the w’s have everywhere positive coefficients. If M,=0, 
the inequalities can be satisfied by taking for the w’s any positive 
quantities. If M, > 0, the first M, — 1 of the w’s can be taken as any 
positive quantities, after which the inequalities impose M, lower lim- 
its on the choice of the last. Since the compatibility of each set of 
inequalities is a sufficient condition for the compatibility of the pre- 
ceding set, the theorem is now established for the case when a and 
6 have no common zero. 

To prove the theorem in general we note first that condition (ii) 
of the theorem is certainly necessary, since if any fiber 7, which is in- 
active for both a and f is removed from the circuit, the removal will 
not affect the activity of other fibers in the system. But this removal 
reduces to a simple chain that part of the J-th circuit which follows 
j,, and its activity is thus uniquely determined by the o’s and hence 
the same for a and £. Therefore, since for such circuits the fiber n, 
must make the same contribution to y, as to z, equations (14) are 
unaffected in this more general case except for the addition on the 
right of these contributions, equations (15) and following are unal- 
tered, and equation (16) still holds when J refers to one of the cir- 
cuits now being considered, provided there is an h, which is a zero of 
one but not of the other of a and 6. Thus the necessity of condition 
(i) of the theorem follows in the general case. 

Suppose now the conditions of the theorem are satisfied and let 
us undertake to determine a o consistent with a and f. First, disre- 
garding the circuits in which a and $ have a common zero, choose the 
o for the remaining system as above. Note that all equations and in- 
equalities in y, z and o are homogeneous, so that the introduction of 
a common positive multiplier throughout will not affect the result. 
Consider any circuit, J , in which a and # are to have a common zero. 
By an interchange of notations, if necessary, we may suppose that 
W=y,,— %, > 0. Then either fiber 1, is active in both a and £, in- 
active in both, or else active in a, inactive in 6. Since the a and 6 
are both admissible, the one case or another can hold only when the 
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threshold of 1, satisfies certain inequalities in terms of the thresholds 
of the other fibers 1x, which inequalities, in fact, limit initially the 
set of AP’s a and f which are to be classed as admissible. But since 
these hold, and since 


Obie Ne hi, am h,,=4, — Zin» 


K referring to a circuit for which the o’s have already been chosen 
and since, finally, the introduction of a positive factor into the y’s, 
z2’s and o’s already chosen will not affect the final result, we may al- 
ways adjust this positive factor so that y,,and z,,have the desired 
signs, and 1, will thus be made active or inactive, as the case may be, 
for each of a and §. But this being done, ao and $, having the pre- 
scribed evenness and oddness, each successive o in the circuit can then 
be chosen in turn to satisfy the requirements as to the activity or in- 
activity of the corresponding fiber, after which, if , is to be made 
active, the value of S, must be increased or decreased, as the case may 
be, in an amount sufficient to cancel the contribution (the same with 
a as with 6) of this fiber. Thus the theorem is proved in complete 
generality. 

The author is indebted to Dr. Alvin M. Weinberg for helpful sug- 
gestions and criticisms on this and the preceding papers in this series. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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The development of a molecular theory of the diffusion drag forces in 
liquids is attempted by considering the liquids as limiting cases of very 
dense gases. An expression is derived for the force on a small particle 
suspended in a nonuniform mixture of such gases on the basis of kinetic 
theory. Another expression for the order of magnitude of the force on a 
larger particle is obtained by introducing certain hydrodynamical consid- 
erations. These results are compared with an expression previously de- 
rived by N. Rashevsky, and estimates are made of the order of magnitude 
of the volume force on the granules and particles within a typical cell 
due to the diffusion of metabolites. It is found that the drag forces, ex- 
erted by a diffusing solute, depend to a very large extent on the physical 
properties of the solvent. 


The drag forces exerted by a diffusing substance must be consid- 
ered in a complete theory of metabolizing systems. It has been sug- 
gested that they may play an important role in cell division and other 
phenomena (Rashevsky, 1938). An expression for the force on a par- 
ticle due to a diffusing substance was derived by N. Rashevsky on the 
basis of asymmetric molecular bombardments. A corresponding ex- 
pression was derived by G. Young (1938) on the basis of a hydro- 
dynamical treatment. In the former treatment the effect of the sol- 
vent was not adequately accounted for, while in the latter, it was neg- 
lected. It is the purpose here to derive a more complete expression on 
the basis of the kinetic theory of gases. 

We shall consider the solvent and solute as non-uniform dense 
gases diffusing through each other in the presence of another kind of 
solute molecules, having only thermal motion. If the stationary solute 
is sufficiently dilute, interactions between these molecules may be neg- 
lected, and we may consider the force on each molecule separately. 
Thus for convenience we shall consider only one such molecule, and, 
in order to distinguish it from the other molecules we shall refer to 
it as the particle. The molecules of the solvent will be referred to as 
class 1, and the parameters for these molecules will bear the subscript 
1. Similarly the solute molecules will be referred to as class 2, while 
the subscript 3 will refer to the particle. 

‘Let N, be the molecular density of molecules 1 of mass m, and 
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diameter o,, for which the most probable thermal speed iS a = 
V2kT/m,. A similar notation will be used for molecules of class 2, 
and for the particle, except that N is not defined for the particle. Let 
c; be the velocity of a particular member of class 7 before, and ¢c; the 
velocity after collision. The molecules of all classes are assumed to 
be elastic spheres. 

Since we intend to extrapolate any results to liquids, we wish to 
consider the effects due to the finite sizes of the molecules. The first 
approximation, f,°, to the distribution function of the velocities of 
class 1 molecules, f; , is given by the Maxwell distribution. An equa- 
tion, whose solution gives a second approximation to f,; = f° + fi* 
for a mixture of two dense gases, is given by S. Chapman and T. G. 
Cowling (1939). Similarly, equations for collision frequency are giv- 
en, so that one may obtain the average change in velocity of the par- 
ticle. The presence of a third type of particle introduces a further 
complication. We shall not attempt this general type of approach, but 
use a procedure in which the effect of molecular size will be considered 
in connection with collision frequency, but neglected in connection 
with the second approximation to the distribution function, f,. The 
effect of the size of the particle on the distribution function will be 
similarly neglected. If only the particle size were to be considered, 
one might treat. its effect on the distribution function as done by P. 
C. Epstein (1924) in obtaining Stoke’s law. The distribution func- 
tions will be considered as independent of position and the effect of 
position will be brought in as is frequently done in this type of prob- 
lem. 

The number of molecules per unit volume having velocities in 
the element, de , of velocity space, is Nfde. Then f,, the distribution 
function for a class 1 molecule, is given by 

(e1-w)2 
fi=a? Us Citta (1) 
where w is the mass velocity of the mixture. A similar expression 
holds for class 2 molecules. For the particle we have 
fpr 2a, 6 (2) 
as we consider the particles to have only thermal motion. 

Consider a collision between a molecule of class 1 and the particle. 
Let k be a unit vector from the center of the particle through the center | 
of the molecule 1. The change in velocity of the particle is given by 
the expression (Chapman and Cowling, 1939, p. 58) 

2m, 


“Sin. ak mace cao (3) 
1 3 
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C31 = C3 — CG, = — Cs. (4) 


If dk represents an element of surface on a sphere of unit radius 
containing the point k , then the probable number of collisions per sec- 
ond of molecules of velocity ¢, with the particle, such that k lies in 
the range dk, is given by 


X31 N, o31? (c;,:k)dk, (5) 
On (o3 ae a1)/2, . (6) 


where y;; is a factor depending on the density of the particles, its val- 
ue being unity for rarified gases (Chapman and Cowling, 1939, p. 
292), and where N, is evaluated at the position of molecule 1 at its 
last collision with any molecule of either class 1 or 2. The net accel- 
eration of the particle is the average change in velocity per unit time, 
or the integral of the expression obtained by multiplying the factors 
from (1), (2), (8) and (5). Multiplying the acceleration by the mass, 
Mz, we obtain F;, as the force on the particle due to bombardments 
of class 1 molecules where 


2m 
F;, = — ee ca I if i N, (es, -k)?fsfskdkde, de; . (7) 
Mx My + Ms Ms 


The value of N, is given by 
Ni =Nu + (onk — Liyets/ Cx) - grad Ni, (8) 


since N, is taken at the point of last collision. The quantity L,. is the 
mean free path of molecule 1 in the mixture of 1 and 2 molecules. We 
shall use the expression corrected for persistence of velocity (Jeans, 
1940, p. 206) in which the y factors have been introduced and in 
which the coefficient a, is taken equal to 0.2 and 6=0.4. Then 


VM, + Me 
Pew O12" Noo X12 VMs sie 62a 01” No Ni VM, + Ms 


The values N,, and N. are the values which N, and N. would take at 
the center of the particle if the latter were absent. Expressions simi- 
lar to (6), (7), (8) and (9) may be obtained for collisions of mole- 
cules of class 2 with the particle by interchanging the subscripts 1 


and 2. 
The factor 7;, depends upon the number densities N; ard Nz. in 


general, and upon o, and «2 as well as upon o;. To a first approxima- 
tion we may write 


Ly = (9) 


st — 3 + a, Ny +-b; No. (10) 
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The value of x; is unity if the effect of the particle is neglected, other- 
wise it is greater than one. The factor 7a: is to be evaluated at do.sk , 
the point between centers at collision. 

Introducing (1), (2) and (8) into (7), and defining v, as a unit 
vector in the direction of the gradient of N,, we obtain for the ex- 
pression under the integral, where o = ai 11: 


N, X31 aes ee 2N, WM, X31 
ieee) e % @ %* (¢e,,-k) *kdkde, de; + aan ae 


C2 C32 


x ff f 6 @ a (¢,- ¥,) (ca k) kd kde, des 


(ys a 5a,N, ia 50,N,) é N, aes aes 
ad a a al SEE I a 031 = y é ete é€ as* 
1 . 


TU Oa Oise 


(11) 


X31 Ly. ON, 


ote Ee Ose 0 Vi 


you mead Ci: 1) 
x if [femens = (cea k) kinden dee 


Cy 


(k 4 v1) (ex . k) 2kdkde, des ee 


In the above expression only terms in the first power of the gradient of 
N, are retained. The mass velocity is considered due only to diffusion 
and hence is proportional to 0N,/dv, = 0N2/d.. For it we shall use 
the expression given by Loeb (1934, p. 263), using somewhat differ- 
ent notations, but in which L,. is given by (9) above: 


i ee (o1 Lag ts 08 Dei) one (12) 
3VaUNe + Ne) OV; 
The second term arises from the first term of the expansion of (1), 
since w, < <a,. The subscript o has been dropped so that y;, , N, and 
N, hereafter are to be taken at the center of the particle. 

From consideration of symmetry, we can see that the net force 
must lie in the direction of the gradient. Thus F,, - v, is the value de- 
sired and the force is then with reference to the direction v,. We thus 
take the component of each term of (11) in the direction v,. The in- 
tegration with respect to dk may be simplified by the substitutions: 


_ b= (k- e533) /¢s; , COS B= (31-01) /C31, M=—k-q, 
where 
q= esc B(v, — COS B es:/€3:), dk =dldm/VW1—m?—F, 
(k-»,) =msinf+lcosp. 
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The integration with respect to / is carried out from 0 to W1-— m? 
as the angle, cos? 1, must lie between 0 and z/2 for a collision to be’ 
possible, while m is integrated from —1 to 1. The integration covers ' 
two of the four quadrants so that because of symmetry the result 
must be doubled. For the integration with respect to the velocity ele- 
ments, a linear transformation of the velocity components may be 
used (Loeb, 1934, p. 551). 

The first integral of (11) vanishes as would be expected, since 
it does not depend upon the gradient, and is essentially the resultant 
force on a particle in a uniform gas without mass motion. The value 
of the integral in the second term of (11) is 


— 2/3 n°? a° a3° Vm, + Ms/V Ms 
so that the complete term becomes 
—— jaa Nevoraitecn, 055 (13) 
The integral in the third term of (11) becomes 
R* O,° ag? (M, + Ms) /3Ms 


so that the complete term becomes 


oN, 


4 
7a aoe kT (y31 = .5a,N, rad .5D1N;) 031° 3 v1 


(14) 
The integral of the last term of (11) was not successfully evalu- 
_ated. Hence two special cases were considered—one in which a, is 
negligible compared with a;, and one in which a; is negligible com- 
pared with a,. In the former case, the integral vanished. In the lat- 
ter, in which the thermal motion is negligible, we obtain for the com- 
plete term 


kTm oN. 
_= * yar Lr 0317 = —, 
0 M4 


15 
We, Se es (15) 


Thus the desired expression must reduce to (15) for large ms but 
must vanish for m; = 0 even when the factor m;/(m, + ms) is re- 
moved. Hence the expression will generally contain m, and m, and 
may be written approximately in the form 


3 N. 3 
— Geet mayen En Oe a fe 
al 3 rl 


where s > l. 
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The net force F;,, in the direction v,, is then given by the sum 
of the terms (13), (14) and (16), or introducing the expressions for 
w, and L,. we have 


__ 160 je Mz (Mm, ets Mz) ok X31 N, O31" 
Hie 27a Mm, + Ms N,+N, 


af 
V 20? N, 41 VM (m, + M2) + 2012” No X12 VM, Mz 


1 oN, 
- y/ Zea N, 2 Yo. VV Mz (m, + mz) + 20127 Nx Yo VM, Me OM (17) 
4 ON, 
— 48 kT (431 + 50, Ny — .50, Ny) 01 oa 
5kT Ms? My + Me Y51 O37 0N, 


3(m, + Ms)* V 20:7 N, Yr. VM + M + 2012? No X12 VD) 0% 


In the present development we are considering the mass velocity 
as due to diffusion only. In general, there may be a mass motion 
which does not vanish with diffusion. If this component is included, 
and the gradients vanish, then the resulting expression gives the force 
on a particle moving relative to the solvent. If m; > > m, and 
o3 > > o» and y;, = 1, this reduces to the expression obtained by P. 
S. Epstein (1924, p. 720) and M. Smoluchowski (1906) for the force 
on a smooth particle which is small compared with the mean free path. 

The force on the particle due to molecules of class 2 in the direc- 
tion of the gradient of Nz is given by an expression similar to (17) 
in which subscripts 1 and 2 are interchanged. The net force due to 
both types of molecules in the direction of the gradient of N. is given 
by F's, — F3,. We shall be especially interested in the case in which 
N, (solvent) is considerably larger than N. (solute) and in which m, 
is considerably larger than m. or m,, where generally m. > m,. Now 
Y12 = Yo. (Chapman and Cowling, 1939, p. 292) and y,, because of the 
above restrictions is approximately equal to y.. The coefficient a, 
gives the amount by which a single class 1 molecule increases collision 
frequency between a class 1 molecule and the particle. This is of the 
order of the volume of the molecule, or because the molecular density 
is of the order of unity, we shall set a, = m,, and similarly a. = Mm.. 
The effect of a class 2 molecule on collisions between class 1 molecules 
and the particle is in general a more complex function, but again the 
effect per molecule, b, , is of the order of ms, at least if the particle is 
not too much larger than the molecules. We shall thus set Ds == ite 
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and similarly b, = m,. Then ys: = ys2. Under these assumptions we 
have for the net force Fs = Fy. — Fi, 


kT 431 o3" ———— a O12 My + mM. 
Fy = — ——*_| (Vm, — 1.12 Vm,) 2-. | 2 
: 24444 N,¢;? Iv : ve) 012° My Me 
(18) 
oN. kT o;? [3 oN. 
ah 16 | a = T 3 - [52 (a5 = o1) a (Mm, — M,) 03 N, | cae 


If, in particular, the molecular diameters of molecules 1 and 2 are 
proportional to the masses, one can show that for m. > m,, and 
a2 > o1, F; is negative, while for m, > m2, 0, > o»,F'; is positive. In 
general, this may not be so. However, from the above equation, one 
would expect a net outward force on particles in a volume element for 
which the divergence of the gradient of solute concentration, m2Nz , is 
negative (the solute being produced in the element) provided the solute 
molecules are larger than the solvent molecules; but one would expect 
an inward force on particles in the element if the solute is being pro- 
duced, if they are smaller than the solvent molecules. This result is 
of interest in the diffusion force theory of cell division (Rashevsky, 
1940). 

The importance of the molecular properties of the molecules in 
determining the force is brought out in (18). It is also clear that the 
effect of the size of the molecules and particle becomes important as 
the density rises, especially since we wish the effect of the difference 
of two net forces. It might be pointed out that the viscosity of carbon 
dioxide was predicted for a range of values including some in the 
liquid state, suggesting that in liquids the momentum transport mech- 
anism is similar to that in gases (Chapman and Cowling, 1939, p. 
209). For a complex liquid like water, such a generalization is more 
doubtful. 

We have considered the molecules as if they were smooth elastic 
spheres. If one considered various types of forces of interaction, dif- 
ferent results would be obtained. In general, the result in such cases 
is very similar to the simplified case. If only numerical factors of the 
order of unity were involved, we would neglect them. But, since we 
are dealing with differences, and the factors might not be the same 
in the two opposing forces, it would be well to investigate under what 
conditions the results would materially differ from the present case. 

Thus far we have had to restrict the diameter of the particle to 
values less than the mean free paths L,. and L.,. If the particle is 
much larger than the mean free paths, we may treat the solvent and 
solute as fluids flowing through each other. G. Young (1940) points 
out that in such a case the force due to the moving fluids is the sum 
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of the forces exerted by each fluid acting independently. If v, and vs 
are the velocities of solvent and solute due to interdiffusion, the par- 
ticle being at rest, the force on the particle due to the solvent and so- 
lute are given by Stoke’s law ; 


F,,=32 1 03 Vi » (19) 
F35 = 3 7 2 O3 Vo. 5 (20) 


The viscosity 2 is taken to be that which the solute would have if it 
were a gas at the same temperature and density. The viscosity 7, has 
a similar meaning, but as we shall consider the case in which 
N, > > Nz, we may also refer to 7, as the solvent viscosity. 

The rate of mass transport, d(m, N,)/dt, across a unit area is 
proportional to the gradient of the mass density and in the opposite 
direction. The proportionality constant is the diffusion coefficient D4. 
of the mixture. But this flow is just the velocity times the mass den- 
sity. Thus 


Nae ’ (21) 
0 Vy 


and a similar expression for the solute. The net force F; is then, as- 
suming N, + Nz is a constant, 


U) U] 
Peewee ie a | oe (22) 
The expression for the diffusion coefficient as given by Einstein holds 
roughly for molecules, thus we may write approximately 


kT 
Dy,.,= ————.. (28) 
32 1 C2 
The expression for viscosity of a gas is 
2y VkTm 
U/ = 3 73/2 o ? (24) 


where y is a factor introduced for dense gases, arising for the same 
reasons as x. The viscosity of a dense gas is only slightly different 
from a rarified gas until the pressure is raised to many atmospheres. 
In the case of carbon dioxide, in which Enskog predicted the viscosity 
into the liquid state (Chapman and Cowling, 1939, p. 289), the range 
of variation is only about five. In the case of water, the range is much 
greater, being of the order of one hundred, indicating that other fac- 
tors may enter, even considering the aggregation of molecules. For 
the sake of analogy, we may write the solvent viscosity in this form. 
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Making use of expressions (23) and (24), we obtain an expression 
for the force in the direction of the gradient of No: 


ee kT a1” C3 “We VMz VM, 0 Ne 
F, = —_——_ ke. (25) 
o2 VM, 


~~, Noon?) N, 0:2 


We are interested in the special case in which N, > > N,. If 
then N,/y, is sufficiently greater than N2/y., we have approximately 
We kT G1" \/Ms 03 ON 


ee — : 
yio® Vm, Ne 02 


(26) 


On the other hand for large enough N., the sign of the force is re- 
versed, and the expression for the force is approximately 


__ kT os ONs 
INS 0. 


(27) 


As we wish to consider the case in which N, > > N,, the solute 
would not be a particularly dense gas, so that continuing to neglect 
the effect of interaction, y. — 1. If the molecular densities are ap- 
proximately the same for solvent and solute, equation (26) becomes 
kTm,1/6 03 ON2 | 


VmMz yi No ys. 


F,= (28) 

Comparing equations (26) and (18) one notes that in (26) the 
force varies only as the linear dimensions of the particle and depends 
on the gradient of the logarithm of the concentration, whereas in (18) 
the force varies with the square of the particle diameter, and with 
the gradient of the concentration. Also the force as given by (26) 
does not change sign when m. < m,, o2. < o,. But from (25), if 
My <<m,, and o. < < o, and Nz is not too much larger than N,, the 
sign of the force may change. 

There are several reasons for the difference in form of the equa- 
tions (18) and (26). In the kinetic theory treatment we consider that 
o; is smaller than the mean free path, while in the second treatment, 
the assumption is made that o; is considerably greater than the mean 
free path so that the fluid may be considered as continuous. This is 
just the reason for the difference in exponent of o, as is discussed by 
P. S. Epstein (1924). The other basic type of difference is that of 
the appearance of the gradient of the logarithm of Nz in (25). But 
it should be noted in the more general expression, from which (18) 
arises, N. does appear in the denominator. Furthermore, if instead of 
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the expression (9) for Li., we use that value obtained by ignoring 
collisions of like molecules (cf. Stefan and Maxwell mass velocity, 
Loeb, 1934, p. 266), a term involving the gradient of the logarithm 
arises. As has been mentioned previously, the term in (18) contain- 
ing the cube of the diameter is doubtful. On the other hand, in the 
derivation of (26) it was assumed that the viscosity, 2, of the solute 
has the value it would have if it were a gas at the same temperature 
and density. 

In N. Rashevsky’s derivation of the diffusion drag force on a 
particle of volume V,, he obtained an expression, 

pa- a V, grad c, 

by considering a continuous pressure function due to molecular bom- 
bardments of a solute of molecular weight, M , and integrating over 
the surface of the particle to obtain the net force (cf. Rashevsky, 
1938, chap. vii). The force on a molecule of the solvent was consid- 
ered in a similar manner. He calculated the non-homogeneity of the 
solvent produced by this force, and showed that this nonuniformity 
results in a force on the particle in the opposite direction to the force 
due to the solute. Thus the above expression was obtained for the 
total force on the particle. However, in considering the solvent mole- 
cules as particles of sufficient size to consider a continuous pressure 
function, the opposing force due to the solvent was underestimated 
and the above a, instead of being about 0.7, becomes practically zero. 
Since some work has been done on the basis of that expression, it 
would be well to compare the order of magnitude of the net volume 
force on many particles which comprise a certain fraction, uw, of the 
solution, for biologically plausible values. In both cases, the net vol- 
ume force is proportional to the gradient. The coefficient of the grad 
¢ from Rashevsky’s equation is of the order of 10’ for u= .1, V, ~ 
a; © 10%, M ~ 10. Using corresponding values in the correspond- 
ing expression from (18), together with values which follow from 
identifying the solvent with water, one finds the same order of mag- 
nitude for the coefficient of the grad c, c = m.N.. The second term 
in (18) is found to be largest for these values. But this term is that 
remaining from the original terms due to the finite size of the par- 
ticle, and is uncertain, since the treatment of the size of the particle 
and molecules was not at all complete. If this term was neglected, 
the resulting coefficiertt of grad c would be of the order of 102 to 104 
instead of 10’, depending on the value of oc; as it varies in the range 


from 10-° to 10-* centimeters, thus materially reducing the magnitude 
of the force. 
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In the theory of forces in cell division, we are particularly inter- 
ested in forces on particles in the range 107 to 10-> cm. We shall esti- 
mate the order of magnitude of the force per cubic centimeter con- 
taining a number of particles which take up a fraction, uw, of the unit 
volume. Taking the fairly plausible values: m, = 6 X 10-2 gm., m: = 
24 10-* gm., 4 —~1 > N; = 107, N. = 10”, (@N-/er2.) = 107, we have 
for the volume force from equation (26) on particles of sizes 10~, 
10-*, and 10° cm. the respective values 108, 10° and 10¢ dynes cm". 
Similarly from equation (18) we obtain the corresponding forces per 
unit volume: 6 X 10%, 3 X 10’, and 3 X 107 dynes cm. The contri- 
butions from the first term of (18) are 10°, 10° and 10‘ respectively 
while those due to the second term of (18) are each 3 X 107 dynes 
em. In view of the previous discussion, the latter term is uncertain. 
In any case, the order of magnitude of the force per unit volume turns 
out to be of the order of magnitude of 10‘ to 10° dynes cm. Using 
the above values for u, M and for the gradient in the expression de- 
rived by N. Rashevsky (1940, chap. iii), one obtains for the force per 
unit volume a value about 10° dynes cm.-?, which is of the same order 
of magnitude or greater than that obtained above. Some of these 
estimates of the volume force are considerably lower than previously 
given, but we shall see in a paper to be published subsequently that 
a lower estimate is still adequate. It will be noted that the discrep- 
ancy between the first and second terms of (18) is smaller when 
forces per unit volume are considered than when the coefficients of 
the concentration gradient are compared. 

Summary. An expression has been derived for the force on a 
particle in a non-uniform mixture of gases. The size of the molecules 
and particle has been taken into account only in an approximate man- 
ner. An attempt is made to extrapolate the results to liquids, and 
using biologically plausible values for the various parameters, the 
force on particles in a unit volume is estimated. A similar expression 
for the force on a particle is obtained on hydrodynamic grounds and 
the latter expression is compared with the former. A similar esti- 
mate of the force per unit volume is made. The two estimates are of 
about the same order of magnitude. These results are also compared 
with the expression obtained by N. Rashevsky. Although his treat- 
ment was not adequate, the order of magnitude of the force per unit 
volume comes out about the same. Thus subsequent work based on 
his equation, especially in the problem of cell division are probably 
not basically affected. 

The author is indebted to Dr. N. Rashevsky for suggesting the 
problem and to Dr. A. S. Householder for checking the calculations. 
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Generalising the previously studied neurobiophysical schemes, con- 
sisting of excitatory and inhibitory elements, a neural mechanism is dis- 
cussed, which may be involved in the perception of combinations of mu- 
sical tones. Equations, giving the total value of central excitation for a 
combination of any two tones are derived, available observations are dis- 
cussed in the light of the suggested theory. 


Since the appearance of the classical work of H. Helmholtz 
(1862), a considerable number of mathematical studies have been 
contributed to the fields of auditory perception (Waetzman, 1912; 
1926). All these studies confined themselves principally to the mechan- 
ism of the peripheral organ. Although even under those restricted 
conditions a remarkable amount of success has been achieved, it is 
doubtful whether a complete theory cannot avoid considering the 
phenomena involved in the central nervous system. The purpose of 
this paper is to make a step in that direction. Although we need not 
commit ourselves to any definite assumption about the functioning of 
the peripheral organs, we shall for convenience put ourselves on the 
point of view of the resonance theory of hearing. This suggests im- 
mediately the picture, that sounds of different pitch are perceived by 
different groups of central neurons, each group excited by one or a 
few adjacent fibers of the basilar membrane. 

The non-linear character of the vibrating structure of the ear 
results in the production of harmonics even when the ear is stimulated 
by a pure tone, therefore any tone results in the excitation of a num- 
ber of distinct neuronic groups, each corresponding to the harmonic 
of the tone. Thus a musical tone results in the excitation of m neuronic 
groups, where 7 is the number of harmonics which correspond to the 
given tone. Theoretically » may be infinite. Practically it is always 
limited, for very high harmonics are of subthreshold intensities. The 
upper limit of pitch which can be heard is also a limiting factor. In 
general, n may be different for different fundamental tones. For sim- 
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plicity we shall consider here the abstract case, where n is the same 
for all tones. In view of what has just been said about the non- 
linearity of the peripheral organ, n is not necessarily determined by 
the tone producing instrument. It might even be completely determined 
by the peripheral organ and the properties of the corresponding part 
of the central nervous system. The number ” may vary with the pitch 
of the tone. Let us, however, first consider the case, that n does not 
depend on the pitch. 

Let two tones of different pitch be presented simultaneously. Let 
the ratio of their fundamental frequencies be p/q, where p and q 
are integers and q<p<2q. The total number of neuronic groups 
corresponding to the two tones is not equal to 27 , because some of the 
groups are common to both tones. The lower tone of frequency » covers 
a range of frequencies », 2y --- iy --- ny. The higher tone covers the 
range (p/q)», 2(p/q)y-:-,1(p/q)y---, n(p/q)y. Of the n harmon- 
ics of the higher tone, those for which i(p/q) > fall outside of the 
range of frequencies of the lower tone. Their number is equal to 
n— (q/p)n=n(p — q)/p. Hence only the number n — n(p — q)/p 
=nq/p of all the harmonics of the higher tone falls within the 
range of the frequencies corresponding to the lower tone. But every 
q-th harmonics of the lower pitched tone coincides with every p-th 
harmonics of the higher one. Hence 1/q of all the harmonics of the 
higher tone, which fall in the range of frequencies » — my, are in- 
cluded amongst the harmonics of the lower. Thus the total number 7’ 
of neuronic centers, corresponding to two simultaneous tones is 


if 2p — 
one Re Se (1) 
qd p Pp 
Denoting 
2p ~ 1 
——————-=27 . 
re (2) 
We find the following values of « for different musical intervals 
Dp Dp 
— x = 
qd q sh 
L 16 5 
1. Minor second iE 1-9375. 4. Major third i 1-80 
: 9 4 
2. Major second 3 1-:8889. 5. Perfect fourth 3 1-75 
Papel 8 roi . 45 
. Minor third P 1-8333. 6. Tritone 39 1-9778. 
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3 16 
7. Perfect fifth 5 16667. 10. Minor seventh *y 1-9375. 
é : 8 : 15 
8. Minor sixth 5 1895 . 11. Major seventh g 1-9333. 
; : 5 7s 
9. Major sixth 3 1-80. . 12. Octave 1 1-5 


If we consider that each neuronic center sends inhibitory fibers 
to all others, then the total net excitation of all the m’ centers may be 
calculated, provided the intensities of excitation of all the groups are 
known. Since at this stage we are discussing the problem purely theo- 
retically, in abstracto, we shall consider here a particularly simple 
case, in which all those intensities are equal. This at first seems to be 
a rather artificial assumption. It may be remarked, however, that 
such a situation will occur if the neurons, leading to the center in 
question reach their limiting excitatory values (Rashevsky, 1938, 
chap. xxil) at rather small intensities of the peripheral stimuli. Un- 
der these conditions, for a very wide range of stimuli, the excitation 
of the n’ centers will remain the same. If the intensity of excitation 
of such a center is J, then since it receives (n’ — 1) inhibitory fibers 
from other centers, its net excitation is: I — b(m’ — 1)I, where b is 
is a constant. Hence the total net excitation is of the form: 


E=A'[n' — b(n’ -1)n'+C'], oS) 


where A’ and C’ are constants, the constant C’ representing a possible 
constant general excitatory (C’ > 0) or inhibitory (C’ < 0): level, 
present at the centers, and A’ is always positive, being proportional to 
I. Because of equations (1) and (2), equation (3) may be written: 


Ha Aton = bier 1) ict Cy. RES) 


The right side of equation (4) considered now as a function of x is a 
parabola having a maximum for a-certain x. Considering that both 
n and n’ = xn are much greater than 1, we may write equation (4) 
as follows: 


E=A(x — Bx? +C) Seater (5) 
where 


, 


A= An: “b= bn: c=f : | be (6) 


For each interval, characterized by given x, we can compute E . Set- 
ting the latter as has been done before (Rashevsky, 1940) propor- 
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tional to the aesthetic satisfaction, derived from the perception of the 
musical interval, we thus may estimate the value of this satisfaction 
and compare it with available rank order measurements. 


FIGURE 1 


Equation (5) shows that F has a maximum for % = &%m = 1/2B . 
Experimentally observed rank order values of the pleasantness of 
simultaneously sounded musical intervals (Valentine, 1913) indicate 
a maximum for the major third and major sixth as shown by the full 
line in Figure 1. To both these intervals corresponds the value 
x = 180. Taking accordingly B = 1/3-6 = 0-2778 , we find in gen- 
eral a distribution of pleasantness values which show the desired 
maxima. The constant A in equation (5) is scale constant, while C 
merely determines the zero line. Neither A nor C affect the rank or- 
ders, which are determined solely by B. With the above mentioned 
value of B we find a general correspondence between theory and ex- 
periment except for the minor and major seconds and the minor and 
major sevenths, for which equation (5) gives too high values. In the 
middle of the musical scale the theory gives too low a value for the 
tritone. The reason for the discrepancy in the case of the tritone 
may be as follows. According to the simple scheme here outlined, a 
mistuning of the interval, no matter how small, results in a tremen- 
dous increase of both » and q and therefore of x bringing the latter 
close to 2 and hence decreasing F. Actually a certain finite threshold 
A of mistuning must be reached before we notice it, although A is 
very small. Neurobiophysically we would expect this, considering 
that two centers, which are too close to each other, may actually over- 
lap and instead of mutually inhibiting they will mutually even en- 
hance each other. For the tritone the value of p/q is very near to 
7/5, the difference between 45/32 and 7/5 being only 0-0062. The 
difference in the ratios p/q for two tones that are one half tone away 
is about 0-1 (v. Hornbostel, 1927). In other words the distance 
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. (82/45) — (7/5) is much smaller than a tenth of a semitone and 
would be very likely to be below the threshold 4. Taking for the tri- 
tone p/q = 7/5, « = 1.8571, we bring the theoretical value to a clos- 
er agreement with observation. A natural reason for the discrepancy 
at the ends of the scale may be seen in the possible variation of the 
mutual inhibition of different centers with the distance between these 
centers (Rashevsky, 1938, 1940; Landahl, 1938; Householder, 1938). 
The geometric distance between the centers corresponding to two har- 
monics will not necessarily be proportional to the distance on the 
scale of frequencies, but is likely to increase with the latter. But the 
frequency distance between the non-overlapping harmonics of two 
tones is the smaller the closer p/q is to either 1 or 2. If the inhibitory 
effect, which is measured by the constant b of equation (3) first rap- 
idly decreases with distance then to remain appreciably constant above 
a certain distance, this will have the same effect as increasing the 
constant B in equation (5) for frequency ratios p/q which are either 
close to 1 or to 2. A very small increase in B gives an appreciable de- 
crease of E.. Making such a correction and putting for the interval 
No. 1-and No. 11, B = 0-2817 and for No. 2 and No. 10, B = 0-2799 
we find a rank order curve as shown by the broken line of Figure 1. 
The full line shows the experimental values obtained by C. W. Valen- 
tine (1913), on 146 persons of both sexes. The rank order curves for 
males and females are slightly different, though show essentially the 
same trend. We find, however, that the calculated values for the uni- 
son and the octave are ridiculously low, way below the seconds and 
sevenths. This may be due to several causes. First of all we made 
for simplicity the assumption of equal excitation of all centers. This 
is very unlikely to happen actually. A center excited by two overlap- 
ping harmonics will be excited much more strongly than a center ex- 
cited by a single harmonic. Since the unison and the octave have the 
most of overlapping harmonics, our assumption will result in too low 
calculated values. 

Furthermore, if we consider the inhibition as decreasing with 
distance, it is to be expected that whereas for intervals which are 
close to the unison and octave, b is going to be larger than for the in- 
tervals further removed, yet for the exact unison and exact octave, 
b will be even slightly smaller than for any other interval. This is to 
to be expected, because the distance between individual harmonics in 
a unison or in an octave, measured on the frequency scale, is larger 
than for any other interval. 

Further improvements of the theory should be made by consid- 
ering the different intensities of excitation for the centers correspond- 
ing to different harmonics. If the law of variation of inhibition with 
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distance is prescribed, the whole problem may be solved graphically, 
formally in the same way as H. Helmholtz (1870, p. 303) constructed 
the effects of beats. Analytically the problem may be handled approx1- 
mately by making b an increasing function of the absolute value of 
the difference y = |(p/q) — (8/2)|, or a function of the linear ex- 
pression a|1 — (p/q)| + B|2 — (p/q)|, where a and £ are constants, 
or some other function of p/q. 

Possible improvements in those directions shall be considered 
in subsequent papers. 

The interest of the above results lies not so much in the any 
agreement of theory and observation which both are still too crude 
as in the circumstance that it shows the applicability to this problem 
of the general neurological scheme, frequently used in mathematical 
biophysics. 

It must be remembered that the rank order of consonance is not 
the same as the rank order of pleasantness. Very tentatively we may 
suggest, that the consonance be measured by the relative amount of 
inhibition at a single center, namely by 1 — bn’ = 1 — bnez. It de- 
creases with increasing «. An alternate possibility however shall be 
treated in separate papers. 
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The standard two factor excitation theories should be called “pre- 
excitation” theories since they apply only to those events occurring just 
up to excitation. A true phenomenological excitation theory which de- 
scribes the whole excitation cycle must involve non-linear equations. The 
nature of these non-linearities is suggested by B. Katz’s subthreshold 
response data. From this data is constructed a “local phenomenological 
characteristic” which is analogous to the current-voltage characteristic 
of a non-linear electrical or mechanical system capable of displaying 
relaxation oscillations. Excitation by constant currents is shown to oc- 
cur where the slope of the characteristic changes sign. The variation 
of the time constant of excitation with degree of response, explained by 
W. A. H. Rushton in terms of a liminal length, is described here in purely 
formal terms. The theory as presented explicitly treats only those events 
in the excitation cycle up to and a little beyond excitation; the complete 
excitation cycle (including recovery and repetition) is mentioned as be- 
ing amenable tc mathematical treatment by an extension of the present 
theory. 


The formal two-factor theory of nerve excitation, (e.g., Rashev- 
sky, 1933) while being adequate to predict time-intensity relations 
for a variety of stimulating currents, is generally assumed to be inad- 
equate to deal with phenomena occurring after “physiological” re- 
sponse has begun. (See, however, Katz, 1936). That this inadequacy 
must be fundamental is apparent from a consideration of the physical 
basis underlying excitation: profound structural changes are known 
to accompany excitation and recovery (Cole and Curtis, 1938), yet 
the excitation equations involve parameters which are constant and, 
therefore, independent of these changes. A possible way out of this 
dilemma has been suggested by several authors (Monnier and Coppée, 
1939; Arvanitaki, 1939; Householder, 1939): viz., that the excitation 
equations are really non-linear, and that the usual linear equations, 
which are valid only during the “physical” response, are first approxi- 
mations to the exact equations. (This statement may be taken to imply 
a definition of physical response: i.e., physical response is that re- 
sponse representable by linear equations.) Such a non-linear theory 
could presumably describe the processes accompanying breakdown and 
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recovery ; however, so far as the writer is aware,! no attempt has been 
made to investigate the specific nature of the required non-linearities. 
In the present paper we shall investigate some possible non-linear 
modifications of the formal excitation theory; the discussion will be 
confined to non-accommodative (single-factor) cases, and for the most 
part, to.sub-threshold effects. 

Phenomenological character of excitation theory. The “state of 
excitation” of a nerve, <, is defined operationally by 


E, — E(t) 
Ey / 


where E, is the resting threshold cathodal impulse (measured either 
in volts or coulombs), and E(t) is the impulse required to just bring 
the nerve to threshold at time t. Before excitation E(t) is the shock 
required to just produce a propagated excitation, and so E(t) is catho- 
dal and e(t) < 1; during the early stages of excitation H(t) is the 
shock required to just quench the excitation, E(t) is anodal, and 
é > 1. Thus, although the intrinsic physical nature of « is a matter 
of conjecture, « can be defined operationally over a considerable por- 
tion of the excitation cycle; over this portion of the cycle it is a per- 
fectly valid variable to use as the basis of a mathematical theory of 
excitation (Monnier, 1934; Schaefer, 1940, p. 166). 

It has been shown (Young, 1941) that if e« is assumed to vary 
according to 


e(t) = (1) 


de 

rT KI =e, (2) 
then « defined by (1) is consistent with (2). Proof of this is required 
because any dynamic equation such as (2) implies a relation between 
e and the current J [or E(t) ] which, in general, will not reduce to (1) 
for impulsive stimuli. If (2) is modified by assuming that k is itself 
a function of e, but K is not, then we have 


d 
Ha kl — kee (3) 


it remains to show that e defined by (3) is consistent with (1). 
Let us define two quantities e’ and e” which vary according to 
de’ 


a ae k's 


1 During the preparation of this manuscript I have had the opportunity ¢ 
: R Y to 
read a very interesting manuscript by Dr. B. Katz in which a se a feel - 
proach to the problems discussed here is given. 1D aD het 
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and 


de" 
—— if os ud ” 
Taga K ke 
where k’ and k” are the maximum and minimum values, respectively, 
of k(e) on the range 0 to e. Then it is clear that 


& Se Se", 


the equality holding, for any arbitrary I(t), only for constant k(e). 
Moreover, for impulsive shocks (Young, 1941), 


eee JT = CG) 
eae E. 


and so e(t) of equation (3) is also given by (1). We have therefore 
proved that the dependent variable of the non-linear excitation equa- 
tion (3) is operationally definable by (1). We shall refer to a theory 
which concerns itself only with variables definable in terms of the 
state of excitation as a phenomenological theory; from this point of 
view, the single-factor theory or its non-linear generalization (3) is 
phenomenological. 

The fact that k is a function of the electrode spacing (Katz, 1939, 
p. 75) would decrease the operational significance of e, unless we 
agree upon some standard electrode spacing; in the remaining dis- 
cussion we shall assume that this standard electrode spacing is main- 
tained. 

Subthreshold effects. Experiments by B. Katz (1937), A. L. 
Hodgkin (1938), R. J. Pumphrey et al (1940), and others on sub- 
liminal responses of nerve have strongly indicated that the linear 
equation (1) should be replaced by the more general equation (3). 
Their evidence may be briefly summarized as follows: after an impul- 
sive conditioning shock, « will attain some value, ¢,, and according to 
(2), will then subside exponentially, i.e., 

de __ may Aes Dae eye 

a He. eee" (4) 
where k is the time constant of excitation. (Note that our & is the re- 
ciprocal of Hill’s k). This statement is the starting point of A. V. 
Hill’s (1936) discussion of excitation. The exponential subsidence of 
the excitatory state is certainly correct if ¢, is small or negative (an- 
odal) ; but, as Katz’s (1937) Figure 4 shows, when «, > .5 marked 
deviations from (4) occur, while for ¢, > 1, ¢ actually increases with 
time. An interpretation of these deviations from linearity in terms 
of a “local” response in accordance with Rushton’s notion of liminal 
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length (1937) has been fairly general (Katz, 1939) ; but from the 
strictly phenomenological viewpoint such considerations are some- 
what irrelevant. All that can be concluded from these curves is that 
the differential equation (4) should be replaced by a non-linear equa- 
tion in which the time constant is itself a function of «, (or, a possi- 
bility which will not be discussed here, K is a function of I). The 
physical interpretation of such a non-linearity is outside the scope of 
a formal theory. For example, from the point of view of a formal non- 
linear theory, it is immaterial whether a variable k is due to an all-or- 
none local response extended over a gradually increasingly area, or a 
graded local response, or (as would appear to be most likely a priori) 
a combination of both. 

From the «(t) curves of Katz’s experiment, the function k(«) 
(which we shall call the instantaneous time constant) can be found 
since, by (4), 


ke) =-= &, (5) 


PROBALY 


& 


FIGURE 1 
The “local phenomenological characteristic”, w(e) vs. ¢, and the instan- 
taneous time constant k(e), from Katz’s data. Indicated points are experimental. 
The dotted portions of these curves are hypothetical. The function k(e) for 
anodal shocks is a constant whose value is uncertain (cf. Katz, 1937); for this 
reason no experimental points are shown on the anodal side. 


and the value of de/dt for each value of « can be ascertained by graphi- 
cal differentiation. If k& is a function of « alone, these e(t) curves 
should consist of a family of “parallel” curves, in which the slope of 
each member of the family is the same for a given value of ¢ . To what 
extent this is true can be judged from Figure 1, where —de/dt vs. € 
and k(e) are plotted from Katz’s data; very much the same results 
are found from the giant axon data of R. J. Pumphrey et al. The scat- 
tering of the points taken from different members of the decay curve 
suggests that k may depend on some variable besides ¢ (e.g., de/dt, e,, 
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or I) or that even in this range K is not constant; however, the scat- 
tering is not bad between e © .5 and e © 1.4 and so we shall suppose 
k to be determined by « alone in the present discussion. 

A semi-analytic approximation to k(e) together with the re- 
sponse curves obtained by using this approximation in integrating 
(5) are shown in Figure 2. A general resemblance to the experimental 
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FIGURE 2 
Local decay curves obtained by integrating (5) with k(e) given by the piece- 
wise linear function on the right. (The positive direction on the k-axis is right 
to left.) The notation h 6, and h #, for the points at which k(e) changes slope is 
suggested by Rushton (1937), Figure 6. 


curves is apparent here; improving the fit by bettering the approxi- 
mation to k(e) is of no particular theoretical value, however. 

The nerve “phenomenological characteristic”. In the non-linear 
theory of electrical oscillations and related phenomena it is customary 
to introduce a “characteristic curve’, a knowledge of which is suffici- 
ent to predict the behavior of the system provided the dynamical equa- 
tions of the system are known. In general, the characteristic gives the 
relation between a current and a potential variable; if the current 
variable is a multiple valued function of the potential variable, the 
characteristic is usually denoted by the letter N (north), while if the 
potential variable is a multiple valued function of the current, the 
characteristic is denoted by S (south) (cf. le Corbellier, 1931, p. 32). 
For systems possessing a linear characteristic, the variational resist- 
ance (defined as the slope of the characteristic) is constant and equal 
to the ordinary resistance. While the system is in a state characterized 
by positive resistance, energy must be supplied from the outside to 
maintain the activity of the system; while the system is in a state 
characterized by negative resistance, the system supplies its own 
energy. From the definition of the characteristic, it follows that the 
variational resistance (or conductance) changes sign at a turning 
point of the characteristic. 
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In attempting to arrive at a corresponding “characteristic” for 
nerve based on phenomenological principles (assuming that one ex- 
ists), we are guided by certain requirements, suggested by the elec- 
trical analogies, which must be satisfied by any curve which plays the 
role of a characteristic in the theory. First, the characteristic must 
be a relation between some variable which behaves like a potential 
(cf. Herrenden-Harker, 1940), and a variable which has the dimen- 
sions of current (in a generalized sense). Second, the slope of the 
characteristic must reduce, in the linear case, to a variational con- 
ductance (or a time constant, in circuits with constant capacitance). 
This means that a theory based on a linear characteristic must be 
identical with the linear theory already known. Third, at the turning 
points of the characteristic, the direction of the energy flow must re- 
verse, and the nerve phenomenon must become self-sustaining. 

We can treat the quantity « as a potential variable without com- 
promising our ideas of potential too seriously. For the equilibrium 
state of the nerve is characterized by the condition « = 0, and the 
value of « is a measure of the deviation from equilibrium. Conse- 
quently, the potential energy of the system (in a strict physical sense) 
can be expanded in powers of «, thus establishing a connection be- 
tween « and the potential energy which is valid even though é¢ is a 
purely phenomenological variable. The quantity « is analogous to the 
displacement variable in a spring near equilibrium. 

On the other hand, the choice of current variable is not so clear. 
The quantity w(e) = kel =e by (5) | satisfies the dimensional 


requirement, and when k is constant, it does give rise to a linear “char- 
acteristic” which yields the original linear theory. However, if we 
plot w(e) vs. « (Figure 1), we find that the turning point of this “local 
phenomenological characteristic” (l.p.c.) occurs almost where the 
physiological local response begins, but certainly not where the exci- 
tation has become self-sustaining. The origin of this difficulty seems 
to lie in the fact that the self-sustaining nerve response is a propagat- 
ed disturbance; the criteria for physiological response and for propa- 
gated response are not identical. This fact was pointed out by W. A. 
H. Rushton in 1937. 

Stability and instability of the characteristic. If we assume that 
K is independent of «, (an assumption which is certainly incorrect 
for e > > 1, and may be incorrect even for smaller €), we can indicate 
the effect of a constant current I by means of a straight line drawn 
across the l.p.c. parallel to the « axis, (Figure 3) ; the line KI intersects 
the characteristic at the points A and B. Point A is seen to be stable, 
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for if the state of the nerve is represented by some point A , and if we 
displace this point toward the right (« increasing), KI — w(e), or by 
equation (8), de/dt, becomes negative and the point returns to A; 
similarly if e is decreased slightly, KI — w(e) > 0 and again the point 
returns to A. On the other hand, a displacement in either direction 
from the point B will tend to increase the displacement (de/dt and 
A « have the same sign) and so B is unstable. The argument present- 
ed here can easily be generalized to show that for steady current stim- 
ulation any point on an ascending branch of the l.p.c. (dw/de > 0) 
is stable while any point on a descending branch (dw/de < 0) is un- 
stable. 

Suppose a graded series of constant stimuli are applied to a rest- 
ing nerve; these may be represented by a family of horizontal lines 
as in Figure 38. From the preceding discussion, it is clear that a con- 
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FIGURE 3 
Stability and instability of the l.p.c. for steady current stimulation. 


stant current, J, will not stimulate if KI intersects the ascending 
branch of the l.p.c. The rheobase will be represented by the line which 
is tangent to the l.p.c., or, from the figure, 

Wo 

= K 
where w, is the maximum value of w(e) and IJ, is the rheobase. From 
a knowledge of the l.p.c. we can predict the response to a constant 
current stimulus: if I > w./K , excitation will occur, otherwise not. 
The response to impulsive stimuli presents certain difficulties if 
we wish to make our interpretation depend on the properties of the 
l.p.c. As we have said, the turning point of the local characteristic 
is known experimentally to occur at about « — .5 instead of e = 1; in 
other words, instability to steady current stimulation occurs before 
instability to impulsive stimulation. Just such a possibility was sug- 
gested by Rushton (cf. his 6, with h6,) in discussing the “liminal 
length”; if stimulation is impulsive, the entire liminal length must be 
excited for instability to set in (i.e., the representative point must 
be brought to the zero of the l.p.c.) while if stimulation is steady, 
instability induced at any one point is sufficient to ensure propaga- 
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tion (ie., the representative point need merely be brought to the 
turning point of the l.p.c.). This picture is very suggestive—but since 
it is based on the assumption that the local response is all-or-none, we 
are not justified in accepting it as the complete explanation of the 
difference in response to impulsive and steady stimuli. 

Thus far we have intentionally neglected any discussion of what 
happens to the representative point an appreciable time after insta- 
bility has set in. To answer this question we would have to know the 
shape of the l.p.c. during the later part of the action. It is difficult to 
assign operational significance to the quantity « during the later stages 
of the action, since, as far as the writer is aware, it is impossible to 
block propagation by impulsive anodal shocks applied at the point of 
stimulation after a certain point in the action has been reached. If it 
eventually proves possible to give an operational definition of e dur- 
ing the subsidence of the action (Monnier, 1934), we would certainly 
find it hard to refrain from extending the characteristic to include a 
second ascending branch (dotted line, Figure 1). Such a character- 
istic (S type), together with suitable assumptions concerning the 
variation of K with e, would enable us to give a phenomenological 
theory of such processes as recovery and repetitive discharge. The 
analogy with systems capable of displaying relaxation oscillations 
would be complete. However, we do not feel justified in going fur- 
ther into these possibilities until some evidence is given for extending 
the characteristic in the manner just described. 

Relations between linear and non-linear theory. It is an empiri- 
cal fact that the linear equation (2) predicts excellent strength-dura- 
tion and voltage-capacity curves. If (2) is only an approximation to 
the non-linear (3), why should this be the case? 

Suppose voltage-capacity or strength-duration experiments are 
performed in which the index of response is not necessarily the pro- 
pagated disturbance ( = 1), but some arbitrary value of the excita- 
tory state, e,, as measured by the test shock method. Voltage-capac- 
ity measurements of almost this type have been performed by H. Ro- 
senberg (1937) but technical difficulties have made the corresponding 
strength-duration experiments unfeasible (although ef. Rushton, 
1932). If es, is small (< .5) so that the linear theory applies, the time 
constant found from such experiments would simply be the initial 
slope of the l.p.c., and would be independent of degree of response. 
As e increases, the linear theory no longer applies: the effective time 
constant will depend on «, and, in general, on the shape of the stimu- 
lating current. 

To determine exactly what intermediate value of & must be used 
in a linear theory strength-duration curve so that it will deviate as 
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little as possible from the non-linear strength-duration curve, we must 
require that the duration of a constant stimulus according to the non- 


linear theory, 
‘| En de 
Wek ie ie(2).é 


and the duration according to the approximating linear theory 
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if Gy ky (en) € 


be as nearly alike as possible, for a given index of response «,. The 
quantities K; and k;(e,) are the approximating linear theory con- 
stants for the index of response «,. An exact solution of this mini- 
mum problem for all currents I(t) seems to be very difficult, and so 
we shall simply write down an approximate solution which is valid 
for large values of the stimulating current: 


Ja (ex) a5 [CeGeds, KK. (6) 


This integral is the first moment of the function k(e). The fact that 
k,(e,) is not independent of J in the next approximation (small val- 
ues of J) may be related to H. A. Blair’s (1936) contention that the 
simple linear theory is inadequate, even for strength-duration data. 
To determine voltage-capacity or strength-duration curves accu- 
rately we should integrate (3) with the function I(t) replaced by the 
particular current used as stimulus. The result of this integration for 
constant currents of various intensities, is given in Figure 4a. (The 
k(e) function of Figure 2 was used here.) These curves are practi- 
cally the same as those given by W. A. H. Rushton (1937, Figure 6). 
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The variation of ¢« under constant current stimulation obtained by integrating 
(3) with the k(e) of Figure 2. The parameters in Figure 2 are arbitrarily taken: 
k, =1,h 0, —.5, K = 1, and the units of t and J are arbitrary. 
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FIGURE 4b 
Strength-duration curves for indices of response ©, — 1, .80, .65, .50. Points 


obtained directly from the curves of Figure 4a are represented by circles. The 
curves were obtained by using the Table 1 values of k,(«,) in the ordinary 
condenser formula (2). 


From the curves of Figure 4a, we can construct strength-duration 
curves for any arbitrary index of response; the curves for «, = 1, 
.8, .65, .5 are given in Figure 4b. The values of k,(e,) determined di- 
rectly from these curves, together with theoretical values of k; (én) 
calculated from (6), are given in the accompanying Table 1. The 
magnitude of the variation in k;(e,) which may be expected is seen 
to be of the order of 100%; this agrees with the results of H. Rosen- 
berg (1937) and of B. Katz (1937). 


Table 1: Values of k,(¢,) from strength-duration curves 
(Figure 4b) compared with values of k,(¢,) from (6). 


en ky, (&,) (Fig. 46) ke, (e,) (Eq. 6) % difference 
3) 1.0 1.0 0 
.65 .86 .92 6+ 
.80 74 79 6+ 

1.0 .b4 .58 T+ 


It is seen from Figure 4b that the non-linear strength-duration 
curves are practically identical with the linear curves which are ob- 
tained by integrating the linear equation (2), after replacing k, by 
k,(e,). This must be considered rather fortuitous, since, by taking 
some different approximating function for k(e), we could get a dif- 
ferent non-linear strength-duration curve. Thus it appears that Rush- 
ton’s prediction of a strict linear theory strength-duration curve, even 
when the liminal length is taken into account, must depend on the 
assumption of an all-or-nothing response at a point. 

Physical interpretations. Although the theory which we have out- 
lined is phenomenological, certain connections with physical interpre- 
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tations seem to be inevitable. As an example, we may ignore the dif- 
ficulties for the moment, and identify « with the condenser voltage, 
V, in a circuit consisting of a condenser, C , shunted by a non-linear 
resistance. For this circuit, assuming no current flow from the outside, 

dV V 

OT AONE 
where 7 is the current through the variable resistance R(i) and, by 
definition, 7R (7) = V(t). Since 7 = ce 
may be determined exactly as we determined w(e)—by graphically 
differentiating a family of voltage decay curves. The V (7) curve may 
be identified with our w(e) provided we plot « vertically and w hori- 
zontally; on the other hand, this characteristic bears little resem- 
blance to the steady state V(i) curve found by K. S. Cole and H. J. 
Curtis (1940) for squid axon. 

The direct measurement of the V (7) characteristic in nerve is 
analogous to the usual direct measurement (with ammeter and volt- 
meter) of the characteristic of a neon tube shunted by a resistance 
and a condenser. However, there is another ‘‘phenomenological” meth- 
od of measuring the characteristic of such a neon tube which, while 
highly impractical, is analogous to the scheme we have used here to 
get at the “characteristic” of nerve (cf. also Blinks and Skow, 1940). 
Suppose, instead of directly measuring the current as a function of 
the steady voltage across the tube, we determine the time course of 
the discharge of the condenser through the neon tube and its shunting 
resistance. This may be done either by actually measuring the time 
decay of the initial charge with a ballistic galvanometer, or, by de- 
termining the extra charge required at each time to just cause the 
neon tube to flash. In this way we can obtain “local excitatory decay” 
curves entirely analogous to those found by Katz, and from these we 
can determine the instantaneous voltage across the tube as a function 
of the current through it. For a simple neon tube circuit this phenom- 
enological method and the previous direct method should give identi- 
cal characteristics. If the experiments of K. S. Cole and H. J. Curtis 
may be considered the proper analogue, in the nerve, to the direct 
measure of the characteristic of the neon tube, then as we have men- 
tioned, the direct and the phenomenological methods in nerve do not 
lead to the same characteristics. The explanation of this difference 
probably lies in the inadequacy of the structural assumptions required 
to make the transition from the phenomenological to the physical 
characteristics. A complete elucidation of this difference will be a ma- 
jor task of any adequate theory of nerve activity. 


, the V(t) characteristic 
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